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ABSTRACT 

Rainfal 1 -runoff data were collected for a 167-2 acre combined 
sewer drainage area at Halifax, Nova Scotia, Thiese data were used to 
compare the performance of the U.K. Road Research Laboratory (RRL) runoff 
model and the RUNOFF block of the U.S. Environmental Protection Agency 
Storm Water Managment Model. 

Both models were found to be capable of producing creditable 
results for a sewerage system that included surcharged pipes and yielded 
runoff from pervious areas in major storms. Calculated peaks were gene- 
rally earlier and higher than observed values. Agreement improved with 
the magnitude of storm events, Hydrographs from the models were more 
similar than were calculated and observed results, especially in larger 
storm events. The EPA RUNOFF block required more computer time than the 
RRL model . 



RESUME 



Les donn^es relatives aux precipitations et au rui ssel lement ont 
Stfe recueillies A Halifax, Nouvel le-Ecosse , pour une superficie de 167-2 
acres drainfee au moyen d'^gouts unitaires. GrSce ^ ces donn^es, on a pu 
comparer les resultats fournis par le moddl du ru i ssel lement R.R.L. 
(U.K. Road Research Laboratory) ^ ceux du modSle E.P.A. (U.S. Environmental 
Protection Agency Storm Water Management Model). 

Ces deux modules ont donnfe des rSsultats plausibles pour un rfeseau 
d'fivacuation fonctionnant en surcharge et ont permis, pendant les gros orages , 
un apport de ru i ssel lement provenant de terrains permeables. Les valeurs 
maximales calcul^es survenaient plus t6t et 6taient plus 61ev§es qu 'en 
r6alit6. L'6cart entre ces deux sortes de valeurs s'est amenuis^ selon 
I 'inverse de 1 'intensity des orages. La similitude entre les hydrographes 
obtenus 5 partir des modules §tait plus grande que celle qui existait entre 
les valeurs calcul^es et les valeurs observSes, surtout durant les plus 
gros orages. La durie du traitement par ordinateur a €t& plus longue pour 
le module E.P.A. que pour le module R.R.L. 
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CONCLUSIONS 

(1) Both the EPA and RRL models demonstrated the capability of producing 
creditable results when they were applied to a system that exhibited 
features they were not designed to represent. 

(2) In most storms, calculated peaks were earlier and larger than 
observed values. Agreement between calculated and observed 
hydrographs improved with the magnitude of storm events. 

(3) Differences between observed and calculated results cannot be 
explained by errors in recorded hydrographs or in subtraction of 
base flows. 

{k) Pervious areas do not contribute to runo^^f, and runoff volume 
varies linearly with rainfall, for most storms on the Halifax 
catchment. In major storms pervious areas make a minor contribution 
to runoff. 

(5) The effect, on hydrograph shape or peak flow, of the inability of 
the RRL method to represent pervious area runoff was obscured by 
the effects of other factors influencing the shape of hydrographs 
calculated in this study. 

(6) Halifax rainfall is not typical of most populated parts of Canada 
and the United States, where rainfall intensities are higher and 
pervious areas may contribute significantly to runoff. 

(7) Use of artificially steepened slopes in the models, for a number 
of pipes known to be surcharged in the actual system, permitted 
both models to function without the effects of surcharging. The 
influence of this assumption on the shapes of calculated hydrographs 
were masked by effects of other factors on hydrograph shapes. 

(8) The resemblance between the results of the two models was more 
pronounced than the resemblance between observed and calculated 
results. The RRL model tended to generate more pronounced peaks, 
with faster rise and recession, than the EPA RUNOFF block. This 
effect was particularly apparent in small storms and in the low 
flow portions of larger events,; 
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(9) Caution should be exercised in assuming values for input variables 
for either the RRL model or the EPA RUNOFF block unless these values 
have been validated by experience on other basins, or by calibration 
of the models for the basin that is under consideration. 

(10) Computer time requirements, and corresponding costs, for running 
the EPA RUNOFF block were considerably greater than the time and 
cost of running the RRL model. Time and cost could be reduced 
by "lumping" of subcatchments. 
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1, INTRODUCTION 

In the period I968 to 1971, the senior author was responsible 
for a study of "Pollution Attributable to Surface Runoff and Overflows 
from Combined Sewerage Systems" (l). That project and subsequent work 
produced storm runoff data for a I67.2 acre combined sewer drainage area 
in Halifax, Nova Scotia. Availability of these data provided an oppor- 
tunity to consider the performance of two runoff quantity models that 
have been developed specifically for urban sewerage systems: the RUNOFF 
block of the U.S. Environmental Protection Agency Storm Water Management 
Model (2), and the RRL model that was developed by the U.K. Road Research 
Laboratory and has been programmed and evaluated by the Illinois State 
Water Survey (3) - 

Computer programs for the EPA and the RRL models were provided, 
on magnetic tape, by Mr. Harry Torno of the U.S. Environmental Protection 
Agency. The EPA program had been adapted for use on the NSTC computer 
facility prior to this study; the RRL program was adapted during this 
project . 



i> DESCRIPTION OF DRAINAGE BASIN 

The location of the drainage basin used in this study is shown 
on Figure 1. The metropolitan region that includes the City of Halifax 
and the adjoining City of Dartmouth has a population of about 250,000. 
The City of Halifax, which occupies 16,000 acres, has a population 
{1971 Census) of 122,000. The peninsula shown on Figure 1, which was 
the City of Halifax prior to annexation of surrounding areas in January, 
1969, occupies 4522 acres and has a population (I97I Census) of 79,255. 
Sewers in the older, peninsular portion of the city are alnxist all combined. 

The drainage basin considered in this study is shown on Figures 
2 and 3. The size of the area is 167.2 acres. The drainage basin 
boundaries are clearly defined, with one exception - an interconnection 
at the junction of Preston and Pepperell Streets in the north east corner 
of the area. No definitive information about the behaviour of flows 
at this junction was available, but observations of flows in several 
storms, and consideration of sewer sizes, slopes and drainage areas, 
suggested that this interconnection would allow a relatively unimportant 
amount of storm runoff to leave the drainage basin at this point. 

Land use in most of the area is residential. Population, 
estimated from the 1971 Census, is 't^?© persons. A commercial area 
(1-3 storey buildings) on both sides of Quinpool Road, occupies 13-5 
acres, 8 percent of the total area. Table 1 shows distribution of types 
of surface in the drainage basin. 

Flow measurements were derived from a critical depth meter in 
a channel installed to measure flows into a retention basin at the outlet 
from the area. Records in 1970 were in the form of a strip chart of 
flow rate. In 1971 this system, which was accurate but unreliable, was 
temporarily replaced by a circular chart recording of flow depth. 
All but one of the storm records used in this study were derived from the 
1971 data. 

Table 2 summarizes meteorological data for the Halifax region. 
Figure k shows rainfall intensity-duration data for comparison with other 
areas. Rainfall data used in this study were recorded, at the location 
shown on Figure 2, by a standard AES tipping buc!<et gauge. 
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FIGURE 2 FAIRFIELD ROAD DRAINAGE BASIN 
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TABLE 1. SURFACE TYPES IN HALIFAX DRAINAGE AREA'- 



% of Area with Each Surface 



Surface 



Commercial 



Res identia 1 



Streets 

Sidewalks 

Paved Walks £ Drives 

Grass Plots (between streets and 
sidewalks) 

Unpaved Walks S Drives 

Lawns 

Roofs (directly connected) 

Garage Roofs 



19.3 
11.2 
16.2 



S.I 
12.2 



15.0 


fs 


IM 


la 


i.f 


^5.0 


20.0 


3.0 



* Residential and commercial land uses. 



TABLE 2. HALIFAX METEOROLOGICAL DATA 
(Long-Term Means) 



Monthly Temperature 

Annual Snowfall (rainfall equivalent) 

Annual Rainfall 

Annual Precipitation 



kk 


.3°F 


7 


9 


inches 


ke 
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i nches 


5't.2 


inches 



10.0 



.1=; 



>- 



a 






1.0 



0.1 



1 — r 




J L_L 



I 



1 



^ SOUTH CENTRAL U.S.A. 

'CENTRAL U.SA. 
^TORONTO 
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^ HALIFAX 
^UNITED KINGDOM 
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FIGURE 4 INTENSITY-DURATION CURVES* 



U.S. curves based on Yarnel 1 data as formulated by Steel (4). 

U.K. curve (MacLean) as given by Watkins (5). 

Canadian data from Bruce (fi) , except Halifax data, which is from 
MacLaren (7) . 



J, RUNOFF MODELS 

The runoff models considered in this study were used, except 
for minor modifications noted below, in the forms in which they were 
received on magnetic tape. These have been described by others (2, 3). 

The programs for both the EPA and the RRL models were developed 
for use on an IBM 360 Computer. Both models have been adapted for use 
on the COC 6'400 machine used in this study. 

Both the RRL model and the RUNOFF block of the EPA model include 
plotted output of computed hydrographs. The RRL model provides a Cal- 
comp plotter output of the computed and observed hydrograph for each 
storm. The EPA output is a line-printer plot of the calculated hydrograph. 
Differences in plotter characteristics made it desirable to write a new 
plotter sub-routine for the RRL program. A plotter sub-routine was also 
written for the EPA program to allow a computer plot of observed and 
calculated hydrographs on a scale that would be the same as the plotted 
output from the RRL model. Figures 5 and 6 are examples of outputs from 
these plotting routines. 

Both models provide for specification of impervious area on the 
individual cards used to describe each sub-area. To save time and 
expense in calibration of the models, each program was modified by provid- 
ing for a C-factor correction, which allowed values on all cards to be 
multiplied by a common correction factor. 

The RRL program provides for a number of input variables (such 
as "street length", "soil group", and "antecedent moisture condition") 
that are not used in the program. Where these were identified they were 
ignored when input data were punched. Table 3 lists values of input 
variables that were assumed in this study. The pipe and drainage area 
configurations used for each model are shown In Figure 7. 

Only storm runoff was simulated by the models. Base flow was 
subtracted from observed results before they were compared with calculated 
hydrographs. 
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TABLE 3. ASSUMED VALUES OF INPUT VARIABLES 

EPA MODEL 

y Percentage of total area that is impervious - ^li% 
m Percentage of impervious area with zero depression 

#* l^°[^^^ - 25% 

if Resistance factors for overland flow- impervious area - 0.013 

.4 - pervious area - 0.250 

Mf Depression storage- impervious areas - 0.039 rn 

- pervious areas - o.l84 in 

:mi Infiltration capacity curve-max. rate - 3.00 in/hr 

-min. rate - 1.00 in/hr 

-decay rate - 0.001)5 sec- I 

RRL MODEL 

1) Percentage of total area that is impervious - 3!*^ 

2) Depression storage on impervious areas - 0.029 in 
3iJ inlet time _ 5'^;^ 
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fts: PROGRAM OPERATION 
k. 1 Runoff Volumes 

Table 3 lists the input variables for which values must be 
assumed in operating each of the models. Assumptions for initial runs 
of the EPA model utilized default values provided by the program. An 
inlet time of five minutes was utilized in the RRL model for all sub- 
basins; this value was derived from calculations for a typical sub-basin, 
based on assumptions used by Stall and Terstriep (3). 

In initial runs of both models, values for directly connected 
paved areas were based on those in Table 1 (81^ impervious in the 
commercial area, kk% impervious in the residential area). Resulting 
values of peak flows and volumes greatly exceeded observed values, and 
the information about contributing areas was re-examined. Of the im- 
pervious surfaces in Table 1, the only ones that can be expected to wholly 
drain directly to sewers are streets and commercial sidewalks. Sidewalks 
in residential areas may drain onto grass plots and lawns, and driveways 
may slope away from the street and drain onto pervious surfaces. Roofs, 
which make up a substantial fraction of the drainage area, are required 
by municipal by-law to be directly connected to sewers. Field visits 
indicated that most roofs were connected, but the amount of water that 
was diverted to pervious surfaces, because of clogged or defective gutters, 
down spouts or plumbing, could not be determined. For these reasons, 
estimates of tributary paved areas based on Table 1 were expected 
to be high. 

Since it was not possible to accurately estimate the fraction 
of impervious area from the characteristics of the drainage area, the 
models were calibrated by comparison with the observed data. Figure 8 
is a plot of storm rainfall vs. observed direct storm runoff. A least 
squares regression yields the relationship R = O.^k (P-0.029), where 
P is storm rainfall and R is total runoff, both in inches, for all data 
except that for the August 11, 1970 storm, which falls to the right of 
the straight line. The linear relationship suggested that, for all 
storms except that of August 11, pervious areas did not contribute 
significantly to rainfall, i.e. the expression can be explained by assuming 



m 



4r- 



AUG.11,1970 



m. 

m 



S 



m 

m 

■em 




0.5 1.0 

STORM RUNOFF (inches) 



1.5 
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that in all storms 3^% of the total area contributed to runoff and that 
all of the water on the contributing area ran off except 0.029 inches 
retained in depression storage. The August 11 storm was an exceptional 
event: the 20 minute rainfall of 1.9 inches/hour almost equals the 10- 
year 20 minute intensity shown on Figure k. In a storm of this magnitude, 
pervious areas might be expected to contribute to runoff if soil infiltra- 
tion capacities were exceeded. 

The RRL model was calibrated, based on the foregoing reasoning, 
by assuming that contributing paved areas represented 3k% of the total 
area and that depression storage equalled 0.029 inches, i.e. 

R = 0.3'*(P-0.029) 
or R = 0.3^{P-0.01) 

where R and P are runoff and precipitation volumes in inches. 

The EPA model considers a fraction of the area (the fraction 
used was the default value of 0.25) to have no depression storage. Taking 
this into account, the corresponding expression is 

R = 0,3k x 0.25P + 0.3^ (1-0.25) (P-0.039) 
or R = 0.085P + 0.255(P-0.039) = 0.3'*(P-0. 01 ) 

Table k compares runoff volumes calculated from the foregoing relation- 
ships with the observed volumes on which the regression was based. As 
expected, the EPA and RRL volumes were identical to each other and to 
the regression results (subject to rounding errors), except for the 
smallest and the two largest storms. In the smallest storm, where 
precipitation was less than the depression storage allowance, the RRL 
model permitted no runoff, while the EPA model allowed runoff from the 
25% of the area with no depression storage. In the largest two storms 
the EPA model permitted runoff from pervious areas; the resulting volume 
exceeded that calculated by the RRL method, which considered paved areas 
only. 

When the August 11 storm was first run on the EPA model, the 
calculated runoff volume was 1.68 inches, significantly greater than 
the value derived from the linear relationship and larger than the 
observed runoff. Reconciliation of calculated and observed values 
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TABLE k. STORM RUNOFF VOLUMES 
(inches in the drainage area) 



Storm 


Rai n 


RRL and Regres- 






Date 


(inches) 


sion Formula" 


EPA 


Observed 






Volume 


Vo 1 ume 


Vo 1 ume 


June 13 


0.028 





0.002it 


0.0025 


May 26 


0.0^*8 


0.0063 


0.0065 


0.009 


June 1-A 


0.056 


0.009 


0.0088 


0.010 


June 22 


0. 176 


0.050 


0.050 


0. Oit^ 


May 28 


0.112 


0.028 


0.028 


0.027 


June 9-A 


0. 128 


0.03^ 


0.03't 


0.032 


May 30 


0. 136 


0.036 


0.0366 


0.032 


June 9-B 


0. 152 


0.0^(2 


0.0^12 


0.036 


June 1-B 


0. 192 


0.055 


0.055 


0.061 


June 23 


0.31 


0.095 


0. 10 


0.090 


June 25 


0.^*3 


0.136 


0.133 


0.130 


June 3 


0.i*8 


0.153 


0. 1^43 


0. 167 


July lii 


0.708 


0.23 


0.232 


0.206 


July 20 


1.63 


0.5^'t 


D.557 


0.52 


Aug. 11/70 


3.86 


1.30 


1.3'. 


1.^3 



>v R = 0. 34(P-0.01) 

required that one of the following assumptions be altered: size of 
pervious area; depression storage depth; or, infiltration rates. To 
account for the difference in calculated and observed volumes the size 
of the contributing pervious area would have to be reduced by 66^, or 
depression storage on pervious areas increased by 0-37 inches. Examination 
of the topography of the drainage area indicated that, in yards and lawns, 
many depressed areas existed from which water was removed only by infiltra- 
tion or evaporation, and many lawns sloped downward away from streets, 
i.e. no direct runoff occurred from a portion of the pervious areas. 
These areas were unlikely to equal more than about 2QZ of the pervious 
area. 

The infiltration rates used in the initial run were based 
on default values provided in the EPA model. The infiltration rate declines 
rapidly with time from the beginning of a storm, reaching a virtually 
constant rate after one hour. In the major storms considered here, 
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rainfall rates did not exceed infiltration rates until two to three hours 
from the beginning of rainfall, when the calculated infiltration rate 
equalled 0.52 inches per hour. It was assumed that the magnitude of this 
value accounted for the high calculated runoff volume. Selection of 
a higher values is supported by references (2) and (9), which cite 
minimum infiltration rates of 0.8 and 1.0 in/hr, respectively, under 
favourable conditions, and by Kelling (10), who measured minitnum rates of 
1.0 in/hr on urban lawns and gardens. The precise effect on runoff 
volume of changing this rate could not be calculated in advance because 
reducing runoff rates in the model reduced detention storage and altered 
infiltration opportunity. Adoption of a value of 1.0 in/hr (Table 3), 
reduced the calculated runoff volume to ^ .3^ inches. Time did not permit 
additional computer runs to refine this estimate. 

Ifri Surcharging 

When the larger storms were run on the EPA model, surcharging 
occurred in two groups of pipes: a run of older pipes at the upstream 
end of the area, and several pipes close to the downstream end. The 
latter pipes produced an obvious effect on the hydrograph in the EPA 
model, which represents the discharge from a surcharged pipe as its 
maximum capacity without surcharging; the effect was to chop the top 
from the hydrograph, the "lost" flow being stored and discharged later 
when the inflow recedes below pipe capacity. To avoid this problem, 
the slope of the affected pipes was increased to ensure that their 
capacity was adequate. Table 5 indicates the changes that were assumed 
for this purpose. In making these changes it was recognized that the 
output hydrograph would be modified somewhat, but it was reasoned that 
the effect would be considerably less than the assumption, implied in 
the EPA RUNOFF bloci<, that the pipe could carry no more than its full- 
f low capaci ty . 

The RRL model calculates flows in a surcharged system by 
increasing the size of those pipes in which inflow exceeds pipe capacity. 
The increased pipe slopes used in the EPA model runs were also used for 
the RRL model to make the modelled systems comparable. In a trial run 
the RRL model was run with the original slopes for the July 20 storm. 
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TABLE 5. CHANGED SEWER SLOPES TO AVOID SURCHARGING 



Pipe No. Slope 



EPA RRL Actual Assumed 



^ 1- 'i .0018 .0)00 

5 1-5 .0008 .0100 

9 1-9 .003^* .0100 

10 1-10 .0034 .0100 

11 l-Il .0033 .0100 

29 1-28 .0008 .0500 
47 2- 1 .0001 .0100 
77 8- 1 .0160 .0300 

30 1-29 .0092 .0400 
tl 1-30 .0139 .0400 



The output hydrograph was identical with that for the system with 
assumed steeper slopes. 

4.3 EPA Model Integration Period 

The EPA RUNOFF block routes surface runoff and pipe flows by 
successive approximations, with the expectation that the result will 
converge to a preset tolerance level in a fixed number of iterations. 
If the calculation fails to converge, this fact will be indicated in the 
program printout. 

Failure of calculated flows to converge occurred in a number 
of pipes as a result of rapid changes in rates of flow in the larger 
storms. Reduction of the integration period used In the calculation 
from five to two minutes in the larger storms (at the expense of increased 
computer time) permitted flows to converge to the established criteria 
in all but a few cases. In the remaining few pipes, examination of the 
printouts indicated that convergence was sufficiently close that no 
significant error would result. 

#,i RRL Model Time Interval 



Initial runs of the RRL model utilized a calculation time 
increment of five minutes. All storms were re-run with a one minute incre- 
ment except the August II storm, for which a two minute increment was used. 
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I, RESULTS 

Figures 9 through 23 present the hydrographs calculated by both 
the EPA and RRL models, together with the observed hydrograph for each 
storm, arranged in order of storm peak magnitudes. As noted above, base 
flow has been subtracted from the observed data to provide a hydrograph 
of direct runoff. Figures 9 through 23 have been traced from the plotter 
output in order to present all of the information in summary form on a 
single figure for each storm. 

Table 6 compares magnitudes and times of occurrence of observed 
and calculated peak flows, and Figure 2k shows observed and calculated 
peak flows. 
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TABLE 6. PEAK FLOWS 
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EPA 






RRL 




Flow Rate 


Time 


Flow Rate 


Time 


Flow Rat 
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Date 




(cfs) 


(hr:min) 


(cfs) 


(hr :min) 


(cfs) 


(hr:min) 


June 


13 


0.8 


0:21 


0.7 


0:24 


_._ 


■— ^.'P^ 


May 


26 


2.3 


2:28 


1.1 


2:20 


3.6 


2:20 


June 


1-A 


2.2 


0:52 


3-2 


1:46 


5.4 


1:43 


June 


22 


5.3 


2:kh 


10.1 


1:54 


9.9 


1:55 


May 


28 


5.3 


2:27 


8.1 


2:14 


9.1 


2:14 


June 


9- A 


7.6 


1:13 


13-6 


1:14 


12.2 


1:16 


May 
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0:45 


7-0 


0:48 


9.2 


0:47 


June 


9-B 
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0:45 


13.6 


0:54 


13.3 


0:56 


June 


1-B 
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1:26 
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0:34 


17.4 


0:34 


June 


23 
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0:33 
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0:32 
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0:33 


June 
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3:45 
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3:26 
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3:29 


June 
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2:45 
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1:44 
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1:47 


July 


li* 
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4:00 
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3:52 
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3:55 


July 


20 


93.0 


3:00 


134.0 


2:48 


137.0 


2:50 


Aug. 
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104.0 


4:30 


144.4 


4:26 


128.6 


4:24 
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S* DISCUSSION 

Figures 9 through 2h and Table 6 summarize the results of 
runoff simulations for 15 storms on the 167-2 acre combined sewer catch- 
ment. Since both the EPA RUNOFF block and the RRL model have been 
calibrated to yield runoff volumes that equal observed values in most 
storms, comparisons must relate largely to the time distribution of 
observed and calculated runoff. 

In most storms, both models provided reasonable reproductions 
of hydrograph shapes, the fit improving for storms of larger magnitude. 
Figure 24 summarizes relationships between relative times of occurrence 
and magnitudes of observed and calculated peak flows in 11 storms. Four 
storms are not included because calculated and observed peak flows 
occurred in different parts of these storms. The figure indicates that, 
in a majority of the storms, calculated peaks were larger and occurred 
earlier than observed values. In smaller storms, as indicated in Figure 
25, calculated peaks exceeded observed values by up to 120^. In larger 
storms calculated peaks were within ^0% of observed values. 

Before concluding that reasons for d i screpencies between calcu- 
lated and observed flows originated )n the structure of the two models, 
the precision of the recorded observations with which model results were 
compared was considered. Observed flows used in this study were believed 
to be within at least 10^ of actual flows in the system. Recorded flows 
up to about 10 cfs were verified by comparison with measurements from a 
second critical depth meter which, during low flows, received the total 
flow from the meter at which flows used in this study were recorded, 
it was not possible to verify the precision of recorded values at high 
flows, but experience and experiment with Pa Imer-Bowlus flumes has shown 
measured values to be within 3% of actual values (8). On a number of 
occasions at high flows, chart readings were verified by manual depth 
measurements. During this study, the sensitivity of observed flows to 
differences in individual interpretation of chart readings was tested 
and found to be within acceptable limits. 

In small storms, base flow values used were of the same order 
of magnitude as storm flows, and errors in subtracting base flows might 
have significantly affected observed volumes and peak flows. Base flows 
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FIGURE 25 MAGNITUDE AND TIMING OF PEAK FLOWS 



29 



from the system considered here were measured during dry weather on many 
occasions, and the results were similar to flows at the beginning of the 
storm events described in this study. Furthermore, errors in estimating 
base flow would not account for distinct differences in the shape of 
observed and calculated hydrographs in several of the smallest storms. 

The models were run for storm flows only, while in the actual 
system storm flows were superimposed on base flow, which was subtracted 
from recorded outflow hydrographs before they were compared with calculated 
hydrographs. The routing procedure used in the model produced a flatter 
hydrograph, with later peaks, if base flow was not included with storm 
flow, i.e. if base flow had been routed with storm flow in the models the 
calculated hydrograph peaks in small storms would tend to be further 
accentuated. 

Resemblance between the results of the two models was greater 
for all storms than the resemblance between the model results and the 
observed hydrographs. The RRL model tended to generate more pronounced 
peaks, with faster rise and recession, than the EPA RUNOFF block. This 
effect was particularly apparent in small storms and in the low flow 
portions of the larger events; in the major storms it was almost unde- 
tectable. The difference presumably arises from differences in the manner 
in which the two models generate inlet hydrographs and route them through 
the system. For example, the way in which depression storage is assumed 
to fill may account in part for hydrograph shape. The EPA model yields 
runoff in any storm, from the time the storm commences, by assuming that 
a fraction of the impervious area has no depression storage. The RRL model 
yields no runoff until depression storage is satisfied, which can eliminate 
or reduce hydrograph peaks if peak rainfall occurs early in a small storm, 
and results in no runoff from storms with rainfall smaller than the assumed 
depression storage (e.g. Figure 9). 

The RRL model does not compute runoff from pervious areas. 
In most of the storms considered here this did not matter because runoff 
from pervious surfaces did not occur. In the August II storm, volumes of 
runoff calculated by both the RRL and EPA methods were less than the 
observed value, because the EPA model coefficients had been altered to 
reduce pervious surface runoff, and because the RRL model yielded no 



runoff from pervious surfaces. The effect of the underestimates on 
hydrograph shape or peaks is unclear because both calculated hydrographs 
were steeper, and peaks higher, than observed values. Studies carried 
out by the Illinois State Water Survey (3) demonstrated the need for 
representation of pervious area runoff in the RRL method, and led to 
development of a version of the RRL method, ILLUDAS, with this capabi- 
lity (9). 

To avoid effects of surcharging in the EPA model in the larger 
storms, slopes were artificially steepened in the 10 pipes listed in Table 
5 and shown on Figure 7. The assumed slopes did not completely eliminate 
surcharging in the August 11 storm; Figure 23 shows a truncated peak on 
the EPA hydrograph, caused by brief surcharging of a major pipe close 
to the outlet. The artificially steepened slopes, and resulting higher 
calculated velocities, tended to yield calculated hydrographs that had 
earlier and higher peaks than actual flows. Time was not available during 
this study to re-run the models with actual slopes for the storms in 
which surcharging did not occur to determine the significance of this 
effect, but It was judged to be unimportant: only three of the pipes in 
the main downstream reach of the system were affected, and the other 
seven controlled a relatively small portion of the total system. 

In larger storms, when pipes were surcharged, hydraulic gradients 
in the real system exceeded pipe slopes, and if the steepened gradients 
remained below street level, system capacity would exceed that for full 
pipes under open channel conditions. Street and sewer slopes in the 
Halifax catchment are such that, with the degree of surcharging encoun- 
tered In the storms considered here, hydraulic gradients should remain 
below street level. Neither the EPA RUNOFF block nor the RRL model was 
capable of reproducing this effect, and the assumed steepening of pipe 
slopes appeared to be the only realistic way to use the models to describe 
the behaviour of this surcharged system. Whether differences between 
observed and calculated hydrographs were the result of this assumption, 
or of other features of the models, remains for further study. Differences 
In the magnitude of observed and calculated values of water released 
from channel storage during periods of recession during the August 11 storm 
(Figure 23) exceeded the volume of the pipes that could be affected by sur- 
charging, and indicated that other factors may be more responsible. 



From the foregoing discussion it is evident that the real 
system considered in this study differed from the systems the models 
were designed to represent, in that the RRL model does not represent 
runoff from pervious surfaces and neither model represents a surcharged 
sewerage system. Despite these differences, in this study both models 
demonstrated enough flexibility to yield creditable results when they 
were applied to a system that exhibited features they were not designed 
to represent. 

The models reproduced flow variations from a system that has 
features common to many mature residential areas that are served by 
combined sewerage systems, including a high proportion of impervious 
surfaces and surcharged sewers. The system is atypical of systems in 
many other parts of North America in that rainfall intensities at Halifax 
are relatively low and pervious area runoff occurs only in severe storms. 
The storm events simulated in this study ranged from those yielding flows 
of the order of magnitude of dry weather flow, to the five year storm. 
Both models represented both hydrograph shapes and flow magnitudes best 
for the larger events, those on which system capacities would be based 
if the models were used for design. Their ability to represent flows 
from minor storm events, which might be important if the models were 
applied or adapted to represent annual or seasonal storm sewage flows or 
combined sewer overflows, was less remarkable. 

In preliminary runs of both models, calculated runoff volumes 
and peak flows were excessive when values of depression storage, imper- 
vious area and infiltration rates were assumed from other sources or 
basin characteristics, i.e. until they were calibrated with observed data. 
This experience argues for extreme caution in assuming values for input 
variables for either model unless the values have been validated by 
experience on comparable basins, or by calibration of the models for the 
basin that is under consideration. 

An important feature of both the RRL and the EPA models is the 
relatively large size of the computer programs. This is particularly 
important in the case of the EPA model. Computer time requirements for 
both models for the large storm of August 11, 1970 were: 
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EPA RRL 

Centra) Processor Time 88.5 sec 7-5 sec 
Peripheral Processor Time 12.6 sec 3.5 sec 
Input-Output 11-5 sec 1.9 sec 

The larger time requirement of the EPA model reflects the difference 
in the way in which the two models route overload and pipe flows. The 
EPA method requires an iterative calculation, employing the Newton- 
Raphson technique, for all basins and all pipes for each calculation 
interval. The RRL employs a routing method that assumes a relationship 
between storage and outflow and allows a direct solution at each time 
step. 

For each model, computer time and cost, and the cost of data 
acquisition and preparation, can be reduced if only an output hydrograph 
is required, and by "lumping" sub-catchments and representing a smaller 
number of pipe elements. Work completed subsequent to this study (11) 
has demonstrated that reducing the number of sub-catchments represented 
by the EPA model, from the Sk used here to 3, could produce virtually 
identical hydrographs for the storms of June 3, June 23 and July \k, 
provided judgement was exercised in defining the width of the lumped 
catchment. 
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